Electric fields generated by quantized vortices 



B. 



A. S. Rukin and S. I. Shevchenko 
/. Verkin Institute for Low Temperature Physics and Engineering, Kharkov, Ukrain^\ 



o 

(N 



5-H 



i 

C 

O 

o 



> 

oo 

00 

On 
O 



X 



It is shown that in a magnetic field quantized vortices in a superfluid obtain a real quantized 
electric charge concentrated in the vortex core. This charge is compensated by an opposite surface 
charge located at a macroscopic distance from the vortex axis. It is determined that the polarization 
caused by the vortex velocity field does not give rise to electric fields outside an infinite cylinder. 
Observation of electric fields created by the vortices is possible only near the end surfaces of the 
cylinder which must be closed with dielectric covers to prevent superfluid leaking. Influence of cover 
properties on the potential created by the vortex is researched. Potential created by the vortices on 
point and ring electrodes are calculated. 



PACS numbers: 67.90.+Z, 67.25.D- 

A series of recent experiments @ revealed that 
flow in superfluid 4 He is accompanied by appearance of 
electric fields in it. Significant efforts [|[ - [13| were un- 
dertaken to understand the nature of the observed phe- 
nomena, but they are unsuccessful up to date. However, 
several interesting results were obtained which concern 
the mechanisms of polarization of both superfluid and 
normal systems. For example, in Mclnikovskiy's arti- 
cle [4j it is shown that accelerated motion of a dielectric 
medium leads to its polarization proportional to the ac- 
celeration. Natsik Q - Q has studied the peculiarities of 
polarization of superfluid systems and, developing Mel- 
nikovskiy's ideas, has found that vortex motion of atoms 
in a superfluid must lead to their polarization caused by 
centrifugal force acting on them. In this case the polar- 
ization vector is directed normally to the velocity of the 
fluid moving around the vortex line, so the polarization 
vectors form a "hedgehog". Unfortunately, observation 
of electric fields caused by this "hedgehog" is rather dif- 
ficult due to their rapid decrease with distance from the 
vortex line (reverse proportional to the cube of the dis- 
tance). The authors [lj] - [l6| found that in a magnetic 
field the vortex line acquires a real electric charge whose 
magnitude is proportional to the vortex circulation and 
is quantized in the same way as the circulation. The 
compensating charge of the opposite sign appears on the 
surface of the system. Macroscopic spatial separation of 
the vortex charge and the compensating one allows to 
measure the electric field caused by these charges and to 
observe the motion of vortices in the superfluid. In this 
message we discuss the significant aspects of measuring 
the electric fields generated by the vortices. 

First we recall the arguments of the articles [lij _ [3 • 
In a magnetic field H a rarefied medium moving with 
velocity v acquires polarization 



v x H 



na- 



(1) 



Here a is the atom polarizability, n is the medium den- 



sity, c is the velocity of light. It follows from the com- 
mon expression for the electric induction of the moving 
medium 



D = E + 4vrP = eE - 



1 



-v x H, 



(2) 



where [i is the magnetic permeability, that the expression 
(|TJ) for the dipole moment P is valid when /i = 1 and 
e = 1 + Anna. 

In a superfluid there exists a "characteristic configura- 
tion" of the velocity field v caused by a vortex, 



M 



V0. 
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Here <j> is the phase of the order parameter. In the case 
of a rectilinear vortex in a uniform magnetic field this 
velocity field generates a polarization charge 



Ppol 



= -divP = 



-n s iJ 2 (rotv s 
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■—H z —s5 2 (r-r v ), (4) 



where s = ±1, ^(r — r v ) is the 2D 6- function in the 
plane perpendicular to the vortex axis. We denote this 
axis by z and assume that the magnetic field is directed 
along this axis. Thus, a polarization charge appears on 
the axis, and its density per unit vortex length equals to 



2-Kh 



(5) 



The compensating charge appears on the surface, and its 
density equals to the normal component of the polariza- 
tion vector P 



Mr) x H] 

a s = an s 

c 

The total surface charge equals to 
an s h 



a s dS - 



Mc 



J (V0 x H) ■ dS. 
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Now we confine ourselves to the case when the vessel with 
superfluid 4 He has the form of a cylinder and the vortices 
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are rectilinear and directed parallel to the cylinder axis. 
In this case the surface element dS can be replaced with 
L z d\, where L z is the vessel length and dl is the vector in 
the (x, y) plane. The total surface charge per unit length 
along the z axis equals to 



an s h 
Mc 



J (V0 x H) n • dl r , 



Mc 



V(f>-dl T = —r7—H z • 2ns. (8) 

v Mc w 



We have taken into account that the phase incursion 
around the vortex equals to 2its. Comparing this result 
with the result ((H) shows that the electric neutrality con- 
dition is satisfied and the total polarization charge equals 
to zero. 

Here we must note that the distance between the vor- 
tex and the vessel wall is in general macroscopically large. 
However, observation of the electric fields caused by the 
vortex is a complicated experimental task. The fact is 
that in the case of an infinite cylinder the electric field 
generated by the dipole moment P from ([1]) is nonzero 
only inside the cylinder. Really, in the absence of exter- 
nal charges the electric induction D satisfies the equation 



divD = 0. 



(9) 



On the other hand, taking into account that in a sta- 
tionary case rotE = 0, we obtain for the components of 
rot D 



rot z D = -4tt H z — — + . , 

c \ ox oy J 
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In an infinite cylinder the velocity field of the vortex does 
not depend on z. In this case rot^D = rot y D = 0. The 
z component rot 2 D also equals to zero because in the 
stationary case divv = 0. Therefore in this case the 
vector D is constant which is equal to zero, as far as at 
infinity outside the cylinder the polarization P is zero 
and the electric field E is also zero. As the result, inside 
the infinite cylinder the electric field is 



E = 4ttP, 



(13) 



and outside the cylinder it equals to zero. Observation of 
the electric field is possible only for a finite sized cylinder 
in the vicinity of its end surfaces. 

Before writing down the corresponding expression, we 
remind that in superfluid 4 Hc a phenomenon called "film 
flow" takes place. Superfluid helium filling a laboratory 
glass begins to climb up its walls and flow out over the 
edge. Therefore experiments on the behavior of rectilin- 
ear vortices in a cylindrical vessel must be carried out in 
a closed system, i. e. the end of the cylinder with helium 



must be closed with a cover. We will assume this con- 
dition to be satisfied and consider separately two cases: 
a) the measuring electrode (which detects the potential 
caused by the vortex) is outside the cover and outside 
the cylinder; b) the electrode is inside the cover. We as- 
sume that the dielectric constant of the cover is e and its 
thickness is d. Presence of a dielectric near a system of 
charges leads to change of electric field in the space. We 
will consider the influence of a dielectric cover (plate) on 
the electric field of a point charge. 

The problem of a point charge e near a dielectric plate 
of finite thickness d has an exact solution [l7|. In the 
case a) the potential outside the plate on the other side 
of the charge equals to 



1-/3 2 



1 - /3 2 exp(-2fcd) 



Jo(kp) cxp(—kz)dk. 



(14) 

Here j3 = Jo is the Bessel function. It is possible to 
calculate the integral when z ^> d. In this case due to the 
factor exp(— kz) significant values of k are of the order of 
1/z. This fact allows to use the expansion exp(— 2kd) rj 
1 — 2kd and to obtain after a simple integration 



<t>(p,z) 



ifzd 
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< 1. 



During the calculation it was also assumed 
It follows from (TTSl) that with the assumptions made here 
the electric field behind the plate coincides in the zero 
approximation with the electric field of a point charge in 
vacuum. The correction increases with increasing e. 

In the case b) the potential inside the dielectric plate 
is given by an expression 



<j>{p,z) 



(1 



1 - P 2 exp(-2fcd) 



x {exp(-fcz) + exp [-2k(d + z') + kz]}dk, (16) 

where z' is the distance from the charge to the surface 
of the plate which is in contact with helium. Integrals in 
(fTBf can be calculated when z <C d. In this case values 
k X/z are significant and exp(— 2kd) <C 1. Expanding 
the denominator in this small addition, we obtain after 
integrating 



e + 1 yjz 2 + p 2 



2 -^4U + (z + 2z> + 2dfY 1/2 + 



(e+ir 



2 -^)rU + (z + 2dy 



-1/2 



(17) 



This result corresponds (in the zero approximation) with 
the result which can be obtained for a dielectric semi- 
infinite space using the electric image method (see e. g. 
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[18j): the potential inside the dielectric is the potential 
of the point charge in vacuum multiplied by a constant 
coefficient j4f. 

Expressions (fT5|) . (fl7|) show that calculating the po- 
tential of a system of charges in both cases (neglecting 
the corrections) is reduced to finding the potential of the 
same system of charges in vacuum. If the electrode is 
mounted inside the cover, the result is multiplied by ^rj. 

First we find the potential created by a vortex on the 
end surface of the cylinder with helium, more exactly, at 
z = +0. 



P(p') •(/>-//) 



( P - pry + z>< 



.3/2 



tfp'dz' 



P(p')-(p-p') , 2 , , , 

jj^dpdz. (18) 

(p-pf + Z^ 



Thus, the problem is reduced to finding the potential 
of a vortex in an infinite cylinder. However, instead of 
finding <fi(p, +0) by calculating the integrals in (fTS!) , it 
is much simpler to obtain the electric field created by 
the vortex. Due to ([13]) the field is E = -4ttP. We 
limit ourselves with a circular cylinder of radius R and 
assume that the length of the cylinder L and the radius 
R satisfy the inequality L 3> R. The velocity field of 
a rectilinear vortex is a sum of the velocity field of the 
vortex line with circulation <fv s dl = 2-kk located at the 
point i"o = (ro, 9q) and its "image" with circulation — 2ttk 
located at r' = (R 2 /r ,9o)- 



r - r o 



r - r„ 



(19) 



where eg and eg/ are unit vectors tangent to circles with 
centers on the axis of the vortex line and its image cor- 
respondingly. In (|19l) - (j2"4")l r and ro are 2D vectors. In 
the Cartesian coordinates this expression equals to 



v s (r) — kz x V In 



r-ro 



(20) 



Substituting this result to ([T]), we obtain the polarization 
P and the electric field E inside the cylinder 



E 



4iran s KH z |r - r | 



(21) 



its "image" on the side surface of the cylinder equals to 
zero. This boundary condition is a consequence of the 
fact that outside the cylinder the electric field and its 
potential are zero, and as far as the potential must be 
a continuous function of coordinates (otherwise, infinite 
electric fields would appear at its discontinuity points), 
on the surface of the cylinder also (j> = 0. Assuming = 
at r = R, where R is an arbitrary radius vector in the 
(x, y) plane drawn from the axis to the side surface of the 
cylinder, we find 



c = _ 4ff an.fcg. ]n ]R-r | 



IR 



Air an sK H Zln R 
c r 



The electric potential created in a semi-infinite cylinder 
on the internal surface of its dielectric cover equals to 



47r an s K,H z 

:+ 1 C 



■In 



r | R 



|r- r bl r o 
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This potential depends on the position of the vortex rel- 
ative to the electrode. Although the change of the po- 
tential on changing the vortex position is weak, it can be 
measured. For liquid helium (density n = 2 • 10 22 cm~ 3 , 
polarizability a = 2- 10~ 25 cm 3 ) in a magnetic field 10 5 Gs 
the potential has an order of magnitude of 10 _8 V. 

Measuring the potential with a point electrode is not 
always a simple experimental problem, especially in a ro- 
tating vessel with helium. Therefore we consider also a 
situation that takes place if the electrode is a thin metal- 
lic ring mounted into the cover in the (x, y) plane. We 
assume that the center of the ring is on the cylinder axis 
and its radius is r. The potential at an arbitrary point r 
on the ring can be expressed in terms of the charges on 
the ring q(r), on the side surface of the cylinder Q(R) 
and the vortex charge q v : 



<Xr,0) 



q(r>) 



Q(R) 



Ir-R| 



dr , + ! ^y^ dK+ i ^iL dTv , 

• - T v \ 

(25) 

Here q(r) is the charge induced on the ring electrode by 
the vortex charge and the surface charge of the side sur- 
face <5(R) = P n (R), where P„(R) is the normal compo- 
nent of the polarization vector, q v (r) = — divP(r) is the 
density of the vortex charge. 

Firstly it seems that we must know the charge distri- 
bution q(r) in order to find the potential on the ring. 
But it turns out that in the case of a conducting ring its 
equipotentiality allows to avoid the calculation of q(r). 
Indeed, using polar coordinates while integrating (|2"5j) . 
one can write down this expression as 



Taking into account that E = — and using (|2T|). we 
find 



4ff «n.*g. ]n ]E Z £Dl +Cj 



(22) 



The integration constant C must be chosen to satisfy the 
condition: the total potential created by the vortex and 



q(0')d6' 



y/2 [1 - cos(6> - 0')} 

Q(6')Rd9'dZ 



^R 2 + r 2 - 2Rr cos(6>' - 9) + Z 2 
q v dz 



\/r 2 +r 2 — 2tqt cos 9 + z 2 



(26) 
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Now we integrate both sides of this equation by 9 along 
the whole ring. In the right side we replace the variables 
9 — 9' = a, and after this substitution the angular vari- 
ables separate themselves. Taking into account that <j> is 
independent of 9, we obtain 



Q(0')d9'R 



da 

\/2(l — cos a) 
dadz 



\Jr 2 — 2rR cos a + R 2 + z 2 



+ 



d9dz 



^/r 2 — 2rro cos 9 + r 2 , + z 2 



(27) 



As far as the whole ring electrode remains electrically 
neutral, the total charge induced on it is zero, i. e. 
/ q(9')d9' — 0, so, the unknown charge distribution q(9) 
falls out of the problem. 

On the other hand, the condition of electric neutrality 
of helium yields J Q(9)Rd9 = —q v . As the result, after 
integrating by z (from — oo to 0) in the two remaining 
integrals in (l27l) we obtain 



2tt0 



q v 7. l-2^COS#+(^ 2 



In 



1 - 2^ cos 



dG. 



(28) 



In order to calculate the integrals in this expression 



Z7T 

1(a) = Jln(l-2acos9 + a 2 )d9 (29) 



we use the following method. Differentiating both sides 
of the equation by a, we obtain 



dl 



-2COS6 1 + 2a 



da J 1 — 2a cos 9 + a 2 



T d9. 



(30) 



vortices the total potential is a sum of potentials created 
by all vortices. 

The information about the distribution of vortices in 
the vessel with helium will be more complete if we use two 
ring shaped measuring electrodes with radiuses i?i and 
i?2 (R2 > Ri)- Analysis similar to the one performed 
above shows that the potential difference between the 
electrodes equals to 



£>2 = q v 



0, 

MJfc/ro), 
\n{R 2 /R 1 ), 



ro > R2] 
R 2 >r Q > Rr, 

Ri > ro- 



If the vortex is outside both electrodes the potential dif- 
ference between the rings is absent. When the vortex 
passes the external ring (it is assumed that the vortex 
moves from the side surface to the axis), the potential 
difference increases and becomes constant and indepen- 
dent of the vortex position when the vortex passes inside 
the internal ring electrode. This passage of the vortex 
into the rings is determined the more exactly, the less 
is the distance between the rings. Thus, installing ring 
shaped electrodes with various radiuses, one can regis- 
ter the quantity, positions and motion of vortices in the 
vessel. 

It is useful to write down the expression for the aver- 
age potential </>(r) which is measured by a ring electrode 
for a system of vortices with identical circulations (and 
therefore with identical charges q v ) , distributed homoge- 
neously in the vessel. The potential in this case equals 
to 

(35) 

Here Q(x) is the Heaviside step function: Q(x < 0) = 
and Q(x > 0) = 1. Introducing the density of vortices n v 
and replacing the summation by with 2D integration, 
we obtain 



Now it is easy to perform the integration by 9. 

dl 7T / . 9 .N 

— = - (1 - asign 1 - a 2 )) . 
da a x ' 

Integrating this equation by a from to a, we find 



7(a) 



0, a < 1; 
2 In a, a > 1. 



(31) 



(32) 



Returning to (|2"51) and taking into account this result, we 
find the potential of the ring with radius r 



In — , rn > r: 
1 In r o < r. 



(33) 



We see that the potential (f> depends on the vortex coor- 
dinate r only if the vortex is outside the measuring ring 
electrode. If the vortex is inside the ring, the potential 
depends only on the electrode radius r. For a system of 



4>(r) = 
= 2nq v n v 



In — I r\dr\ 
r 



Lo 



In — I r\dr\ 
ri 



Trq v n v (R 2 — r 2 ). (36) 



It is obvious that the potential reaches its maximum 
value in the center of the vessel, i. e. at r — 0. Homoge- 
neously distributed vortices form a lattice (most probably 
a triangular one). The potential </> obtained above is the 
value of potential averaged by all possible positions of a 
lattice with given density n v . 

As it is well known, rectilinear vortices are formed in a 
cylindrical vessel when it rotates with angular velocity f2 
greater than the critical velocity fi c = (h/MR 2 ) ln(i?/£), 
where £ is the coherence length. If 57 ^ fl c , the quantity 
of vortex lines will be very large, and in the 57 — ¥ 00 limit 
they imitate rotation of the superfluid component as a 
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whole with angular velocity Q. The condition rotv g 
gives the density of the vortices 

n v = —-. 

irh 



2Q 



(37) 



Substituting this expression for n v into Q36p and replacing 
q v with its value from ([5]), we obtain the result 



(r) = 2 -^^n(R 2 



(38) 



We can get exactly the same result if we assume that a 
liquid with density p s performs solid rotation with veloc- 
ity v s = f2 x r, find the polarization ([1]) for this velocity 
field, the corresponding electric field E = — 47rP and in- 
tegrate this field from r to R. Thus, the theory developed 
above, which operates with separate quantized vortices, 
gives a correct description of the system's behavior in the 
continuous medium limit. 

Until now we considered only the supcrfluid compo- 
nent. But at nonzero temperatures there always exists 
the normal component. In a vessel rotating with angular 
velocity VL the normal component performs solid rotation 
with velocity v = fl x r. Calculating the polarization ([1]) 
for this velocity field, we find the addition to the electric 
field and its potential 4>{r) caused by the normal compo- 
nent. The addition to the potential is obtained from (|3"8")l 
by replacing the superfluid density n s with the normal 
density n n . Therefore, to observe the effects caused by 
separate vortices, one must carry out the experiments at 
temperatures sufficiently low to make the normal density 
much lower than the superfluid density. For He II this 
condition is satisfied below 1 K. 



It is more convenient to measure not the static po- 
tential, but its variation in time. When a vortex passes 
inside (or outside) a pair of measuring ring electrodes, 
the potential difference between them undergoes a jump. 
If more than one vortex passes, the change of the poten- 
tial is correspondingly multiplied by the quantity of the 
vortices. 

Let us briefly consider the behavior of vortices in a 
rotating vessel. Using the known expression for the free 
energy of a vortex in a rotating vessel [l9j] , one can obtain 
that at angular velocity lo > lo\ = nj R 2 the energy has 
a maximum at a distance xq = y' R 2 — k/lo from the 
axis and a minimum on the axis. Vortices formed near 
the walls of the cylinder try to overcome the potential 
barrier and accumulate on its external side. Therefore 
it is reasonable to place the measuring electrodes with 
radius Xq. 

To observe the motion of vortices in the vessel and not 
their static distribution, one must rotate the vessel with 
variable angular velocity. Numerical calculations show 
that at low angular velocities one to six vortices locate 
symmetrically around the cylinder axis, i. e. on a circle 
with its center on the axis and radius approximately one 
half of the cylinder's radius. When the angular veloc- 
ity (and therefore the quantity of vortices) increases, the 
radius of this circle also increases, and one more circle 
appears inside it, then the process is repeated. Placing 
a pair of measuring electrode with appropriate radiuses, 
one can observe large quantities of vortices periodically 
going out of the rings. 

However, these problems require more detailed analysis 
to which we suppose to revert later. 
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